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ABSTRACT 



Context. The nearly edge-on galaxy NGC 4945 is one of the closest galaxies where an AGN and starburst coexist, and is one of the 
brightest sources at 100 keV. Near and mid-infrared spectroscopy have shown very strong obscuration of its central region, rivaled 
only in strength by some of the most deeply obscured ULIRGs. 

Aims. Determine the spatial distribution of ISM emission features in the central 426x426 pc^ of NGC 4945. 

Methods. We map the central region of NGC 4945 in three of the four Spitzer-IRS modules (SH, SL and LL). In particular, we 
produce maps of the flux distribution of the starburst tracers [Ne II], [Ne III], [S III] and [S IV] at 12.81, 15.56, 18.71 and 10.51 //m, 
respectively, and a map of the AGN narrow-line region tracer [Ne V] at 14.32 /jm. In addition, we mapped the S(l), S(2) and S(3) 
pure rotational lines of H2, which trace the distribution of warm molecular hydrogen. Finally, we obtained an extinction map (Ay) 
based on the apparent strength of the 9.7 /jm silicate absorption feature. 

Results. At a spatial resolution of ~5" our extinction map traces the contours of the starburst ring. The highest extinction is, however, 
found at the nucleus, where we measure Av(9.85 ;um)«60. Within the uncertainty of the astrometry all emission lines are found to 
peak on the nucleus, except for the warm molecular hydrogen emission which shows a maximum 60-100 pc NW of the nucleus. 
We favour a scenario in which the lower H2 0-0 S(l) and S(2) rotational lines originate mainly from an unobscured extra-nuclear 
component associated with the super- wind cone observed in the HST NICMOS map of the H2 1-0 S(l) vibrational line. For the 
[Ne V] emission we infer an attenuation of a factor 12-160 (Ay =55- 112) based on a comparison of the ratio of our [Ne V] flux and the 
absorption-corrected 14-195 keV Swift-BAT flux to the average [Ne V]/BAT ratio for Seyfert 1 nuclei. The high attenuation indicates 
that [Ne V] and [O IV] cannot be used as extinction-free tracers of AGN power in galaxies with deeply buried nuclei. 

Key words. Galaxies: active - Galaxies: nuclei - Galaxies: individual: NGC 4945 - ISM: lines and bands 



1. Introduction 

At a distance of ~3.82 IVIpc (~18.5 pc/arcsec, 
iKarachentsev et al.1 l2007h the active galaxy NGC 4945 is 
one of the closest galaxies that host both, an AGN and starburst. 
Earher X-ray observations showed ev idence for a hidden AGN 
dlwasawa et al.ll993tlGuainazzi et all2000.) . These observations 
revealed a Compton-thick spectru m with an absorbing c olumn 
density of Nh = 2.4 X lO^^ cm'^ (iGuainazzi et"ani2000i) . The 
nucleus of NGC 4 945 is one of t he brightest extragalactic 
sources at 100 keV (iDone et al.lfr 996), and the brightest Seyfert 
2 AGN at > 20 keV ft toh et a'il2008.) . 

Rather than a point source marking the presence of the AGN 
in optical and near-infrared images, the nearly edge-on (/ ~ 80°) 
line of sight to the central region reveals evidence for strong and 
patchy extin ction, which is especi ally apparent in HST-NICMOS 
H-K maps ( iMarconi et al.ll2000h . A dust lane aligned along the 
major axis of the galactic disk obscures parts of the central re- 
gion just southeast of the /T-band peak. The /T-band peak itself 



* We dedicate this paper to the memory of our esteemed colleague 
and advisor, Alan Moorwood (1945-2011), who pioneered the infrared 
spectroscopic study of NGC 4945. This work is based on observations 
obtained with the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under NASA 
contract 1407. 



lies ~ I" west (bu t, within the uncertain ty) of the position of the 
H2O mega maser (i Greenhill et al.|[l99 7). which we adopt as the 
location of the AGN. The mega maser emission allows to esti- 
mate the mass of the central black hole, Mbh ~ 1-4 x 10^ M©, 
which together with its inferred 0.1-20 keV bolometric lum i- 
nosityof~2xlO'^3ergs ' (~ 5x10*^ L^. IOuainazzi et al](l2000h ) 
indicates that the central engine radiates at ~10% of its Edington 
luminosity (Gr eenhill et al. 1997; Madejski et al. 2000). 

Estimates from IRAS observations indicate that about 75% 
of the total infrared luminosity of the galaxy (Ljr = 2.4 x 
10'" Lq) is generated within an elongated region of 12" x 9" 
(about 222 x 167 pc^) centered on the nucleus dBrock et alj 
19881) . The structure of this region, as shown in high resolu- 
tion HST-NICMOS observations of the Paa line, is consistent 
with a nea rly edge-on starburst disk with a 4.5" (~ 83 pc) radius 
(iMarconi et al...2000.) . embe dded in a r == 8" ( ~150 pc) inclined 
(62 degrees) molecular disk (IChou et al.ll2007h . 

Although the star formation and supernova rates in the nu- 
clear region were originally estimate d to be moderate (~ 0.4 Mp 
yr ' and ~ 0.05 yr ', respectivelv: iMoorwood & 01ivalll994l) . 
more recent estimates based on high resolution (angular scale of 
0.3 pc) radio observations and estimates of supernova remnant 
source counts, sizes and expansion rates, lead to a type II super- 
nova rate of > 0.1(v/10'*) yr"', and star formation rate (SFR) 
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limits of 2.4(v/10'*) < SFR(M > 5 Mq) < 370 Mq yr"', where 
V is th e shell radial expansion velocity in km s"^ (iLenc & Tin^avl 
l2009i) . These supernova and star formation rates are, within 
a fact or of two, similar to those estimated in NGC 253 and 
M82 (iPedlar et al.ll2003t iLenc & Tingavll2006h . which are also 
nearly edge-on starburst galaxies, have similar dista nces close 
to 4 Mpc and all show similar infrared luminosities dRice et alj 
11988) . The impact of the central starburst on the circunmuclear 
region is large, as revealed by the pre sence of a conical cavity 
evacuated by a supernova driven wind dMoorwood et al.iri996bl; 
[Marconi et al. 2000). The HST-NICMOS images of the Hj 1-0 
S(l) line at 2.12jum show that the edges of the cavit y extend out 
to 5" (~93pc) north from the NGC 4945 nucleus (.Marconi et alj 
120001) . 

Near and mid-infrared spectroscopy indicate that the ISM 
line of sight to the central region of NGC 4945 is very different 
than that of other nearby starburst galaxies, such as M 82 and 
NGC 253, which present similar inclination angles. The strong 
absorption features of both volatile (CO and CO2) and refrac- 
tory (H2O) ices observed in the 2.4-5 //m ISO-PHT-S spectrum 
indicate the presence of shie lded cold molecular clouds obscur- 
ing the NGC 4945 nucleus dSpoon et al. 2000, |2003|) . The de- 
tection of a strong absorption feature of XCN ice at 4.62 jum by 
I Spoon et al] feOOO, 2003) indicates that these molec ular clouds 
have been processed in an energetic environment dLacv et alj 
11984) . and processed ice is suggested to be a common charac- 
teris tic of dense molec ular material in star forming galactic nu- 
clei dSpoon et alj|2003l) . A highly distorted PAH emission spec- 
trum produced by a very deep 9.7 jum silicate absorption fea- 
ture is more evidence f or the unusually strongly obscured nu- 
clear region dSpoon et al...2000: Brandl et al...2006.) . This makes 
NGC 4945 a unique nearby laboratory to study an environment 
that can be found only in distant ULIRGs. 

Previous mid-infrared spectroscopic observations could not 
confirm the presence of an AGN in NGC 4945. Just an up- 
per limit of the 14.32 /^m [Ne V] emission, considered a 
tracer of AGN narr ow-line regions (e.g.. lMoorwood et al.ll996at 
[Genzel et al."l998'), was obtained from ISO-SWS observations 
(Spoon etal. 2000). The VLT-ISAAC observation of the 3.93 
//m [Si IX] line, commonly observe d in the soft X-ray photoion- 
izedgas of many Seyfert galaxies ( Oliva et al. I 994 : iLutz et all 
I2OO2I) . resulted in a non-detection dSpoon et alj|2003h . Only re- 
cent observations with the more sensitive IRS-SH spectrograph 
on Spitzer allowed the detection (although at a very faint level) 
of the 14.32 ^m [Ne V] l ine towards the NGC 4945 nucleus 
dBernard-Salas et alJl2009l) . The [Ne V]14.3yum/[Ne II]12.8jum 
flux ratio found in NGC 4945 is ~ 0.007, which indicates an 
A GN contribution of l ess than a few percent using the diagram 
bv lFarrahetal] d2007h (their Fig. 16). This detection of [Ne V] 
relative to [Ne n] is about 10 times weaker than what was ob- 
served in other Seyfert galaxies like Mrk 266 and NGC 1365 
(Hirnard-Salas et al. 2009). 

In this paper we study the mid-IR properties of the nuclear 
region of NGC 4945. Spitzer-IRS spectral mapping observations 
of a 23" X 23" (about 426x426 pc^) region are presented. The 
mapping capabilities of Spitzer-IRS allow the study of a num- 
ber of properties of the nuclear region in NGC 4945. For in- 
stance, the extent of the AGN coronal line region as traced by 
the [Ne V] line, the extent of the crystalline silicate absorbing 
region (and a search for the sources responsible for the presence 
of crystalline silicates), disentangling PAH emission and silicate 
absorption along the line of sight, and the excitation and age of 
the circumnuclear-nuclear starburst based on line ratios of for- 
bidden lines. Although we actually used the Spitzer-IRS mod- 



ules SH, SL and LL, in this work we present the most interesting 
results from the SH and SL spectral maps only, since we focus 
on the analysis of the spatial distribution of ISM emission and 
absorption features. We present for the first time the maps of the 
starburst tracers [Nell] 12.81yum, [Nelll] 15.56 jum, [S III] 18.71 
jum, and [S IV] 10.5 1 jum, and a map of the AGN narrow-line re- 
gion tracer [Ne V] at 14.32 /vm, and analyze the ratios between 
some of these tracers as different diagnostics. The organization 
of this article is as follows. In Sect. [2] we describe the observa- 
tions and the data reduction. The maps obtained are presented 
in Sect. [3] The analysis of the data is presented in Sect. H] The 
conclusions and final remarks are presented in Sect. |5] 

2. Observations and data reduction 

We mapped the cent ral region of NGC 4945 with the InfraRed 
Spectrograplfl (IR S -iHouck et alj|2004 on board of flie Spitzer 
Space Telescope dWerner et al.l 12004*) through the guaranteed 
time observation (GTO) program P40018 (PI: H.W.W. Spoon). 
The spectral maps were done in moderate resolution (R~600) 
with the Short-High (SH; 9.9 - 19.6 yum) IRSQ module, and at 
low resolution (R~60-120) using the Short-Low (SLl; 7.4 - 
14.5 //m, SL2; 5.2 - 7.7 ^m) and Long-Low (LLl; 19.5 - 38.0 
jum, LL2; 14.0- 21.3 //m) IRS modules. 

The SH module was used to measure mainly the fine- 
structure emission lines (starburst tracers) [S IV] at 10.51 jum, 
[Ne II] at 12.81 yum, [Ne III] at 15.56 jum, [S III] at 18.71 jum, 
the AGN narrow-line region tracer [Ne V] at 14.32 jum, and the 
molecular hydrogen (pure rotational) lines, H2 0-0 S(2) and H2 
0-0 S(l) at 12.3 yum and 17.0jum, respectively. With the SL mod- 
ule we obtained maps of the H2 S(3), and the silicate absorption 
feature at 9.7 jum. 



2.1. SHmap 

The IRS-SH map was obtained in July 2007 (at the beginning of 
the summer visibility window) when the SH slit was relatively 
aligned with the minor axis of NGC 4945. The area mapped was 
limited to a set of 3 parallel by 10 perpendicular pointings cen- 
tered on the nucleus, and covering an area of 21" x 26". The top 
panel of Fig. [T] shows the orientation of the SH slit overlaid on 
the Spitzer-MIPS 24 pm map of the central region of NGC 4945 
(by courtesy of Varoujian Gorjian, private communication). 

The quality of the map was optimized by splitting up the 
mapping into 4 1 -cycle maps of 120 sec, each of which was al- 
ternated with a staring mode off-source sky measurement to keep 
track of rogue pixels and to obtain an accurate measure of the sky 
background. The total time used for the SH map, including sky 
positions, was 6.5 hrs. 



2.2. SL map 

Due to a Spitzer anomalyQ that began a few hours before the 
original schedule of the observations, the IRS-SL map of the nu- 
clear region of NGC 4945 was delayed till August 2008, when 
the SL slit was no longer perfectly aligned with the galaxy ma- 
jor axis. The observations were done with 5 cycles of 14 sec 



' The IRS was a collaborative venture between Cornell University 
and Ball Aerospace Corporation funded by NASA through the let 
Propulsion Laboratory and the Ames Research Center. 
' http://ssc.spitzer.caltech.edu/irs/highleveloverview 
^ http://www.spitzer.caltech.edu/news/856-feature06-25-Engineers- 
Studying-Spitzer- Anomaly 
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Fig. 1. Spitzer-MIPS 24 yum map (Varoujian Gorjian, priv.comm.) of 
the central region of NGC 4945. The galaxy major axis runs vertical. 
North is toward the upper right, East toward the upper left. Along the 
declination axis the grid spacing is ~1 arcminute, and about 12 sec- 
onds along the R.A. axis. The region mapped with the IRS-SH module 
(top panel) is about 21" x 26", while a larger region of 37" x 67" was 
mapped with the IRS-SL (bottom panel) module. The red and yellow 
grids correspond to the SLl and SL2 spectral orders, respectively. The 
sky (off-source) positions observations are indicated by the grids to the 
left and right of the galaxy major axis. The images were produced with 
the Spitzer/Leopard software package used to query the archive and 
download the data. 



ramps, which allowed to map a region of about 37" x 67" in 2 
parallel by 19 perpendicular pointings centered on the nucleus 
of NGC 4945. 

In order to sample the background continuum and to identify 
rogue pixels, we also obtained off-source spectra with the same 
integration times and number of cycles at two positions suffi- 
ciently above and below the galaxy disk. The off-source spectra 
were taken before and after the SL mapping. The total time used 
for the SL map, including sky positions, amounts to 3. 1 hrs. The 
orientation of the SL slit, overlaid on the Spitzer-MIPS 24 fim 
map of NGC 4945, is shown in the bottom panel of Fig. [T] 

2.3. Data reduction 

The Basic Calibrated Data (BCD) were pre-processed with the 
Spitzer pipeline versions S16.1 for SH, and S18.1 for SL. Rogue 
pixels were cleaned, and maps were built, using the IRS mapping 
reduc tion package CUB ISIvfl designed by the SINGS0 legacy 
team (ISmith et alJl2007h . With this tool we performed the flux 
calibration, background subtraction, and estimated the statistical 
uncertainty at each spectral wavelength. We create the first spec- 
tral cubes from the slit observations combined with CUBISM. 
After creating the first cubes, we furthe r cleaned the data u sing 
the back-tracking procedure described in lSmith et alj (|2007|) . All 
the pixels with flux uncertainty larger than 50% were flagged, so 
they were not used in subsequent reconstructions of the spectral 
cubes. We iterate on the cleaning procedure by performing visual 
inspections throughout the maps, of small sections of the spectra 



around the emission lines, PAH and silicate features of interest, 
and we reconstructed the cubes after cleaning new pixels. 

The spectral cubes obtained with CUBISM are in units 
of surface brightness, MJy sr ', but we convert them to 
W cm"^ jum"' sr ' to work and present the spectra. We use the 
units of surface brightness W cm"^ sr ' to present the maps, and 
flux density W cm"^ for the integrated total fluxes. For this we 
use the conversion 1 arcsec^ = 2.3504 x 10"" steradian, know- 
ing that the pixel sizqjof the IRS/SH map is ~2.26 arcsecs /pixel, 
which leads to 1 .2005 x 1 "* steradians that we need to multiply 
by in order to get the maps in units of flux density. In the case of 
IRS/SL the pixel size is ~1.85 arcsecs/pixel, so we have slightly 
smaller 8.04424 x 10"" steradians per pixel than in the SH map. 



3. Analysis and results 

During the data reduction and spectral analysis process we found 
that the SH spectral order{] are mismatched at levels that vary 
across the mapped region. This effect is smoothed (although not 
completely solved) by increasing the overlap between the orders 
(in the wavsamp calibration file) and by increasing the size of the 
aperture used to extract 1-D average spectra, as discussed below. 
We also observed module to module mismatches, which inhibit 
us from reliably combining the SH and SL modules to perform 
further analyses like, for instance, using the H2 S(3) line as an 
extinction indicator These mismatches in the SH orders and be- 
tween modules have not been addressed beforql. probably be- 
cause they are apparent only in high S/N data as those presented 
in this work. The IRS-SH/SL orders/modules mismatch problem 
seen in high S/N data is a new issue that we are still studying to- 
gether with members of the SINGS team. Fortunately, most of 
the emission and absorption features we are interested in are 
found within the spectral orders, and those features observed 
near the edge of the orders are treated with extreme care or not 
addressed at all. 

The average spectrum, considering the whole field of view 
(FOV) of the IRS/SH map, and the average spectrum of the 
IRS/SL map (extracted from a similar FOV as that of the SH 
map) are shown in the top and bottom panel of Fig. |2] respec- 
tively. The most important fine-structure emission lines and PAH 
features are indicated in both spectra. 

We developed our own IDL procedures to work on the 1-D 
spectra that were used to build the final maps. These 1-D spectra 
were extracted from the clean SH cubes by averaging the spec- 
tra over a moving 2x2 pixel aperture. The whole field of view 
was covered with a step size of 1 pixel, so the reduced map con- 
tains only 10x10 pixels of ~ 2.3" each (where 1" corresponds 
to a: 18.5 pc at the distance of 3.82 Mpc) but covering the same 
area as the original map. However, the last slit observation of 
the IRS/SH map covers just about one third of the spatial region 
associated with each pixel. Therefore, we conclude that the last 
column of the original 11x11 grid map produced by CUBISM 
is not representative of the actual spatial region associated with 
that slit column, and we do not use it in our maps. So the final 
maps reduced with the 2x2 aperture results in 9x9 pixel maps, 
covering a slightly smaller region than the original maps. The 
coordinates assigned to each pixel of the resulting map are those 
of the center of the 2x2 aperture, and the pixel fluxes reported 



http://ssc.spitzer.caltech.edu/dataanalysistools/tools/cubism/ 



* See http://coolwiki.ipac.caltech.edu/i ndex.php/Unitsl for a descrip- 
tion of Spitzer units and conversions 
' http://ssc.spitzer.caltech.edu/irs/irsinstrumenthandbook/ 



http://ssc.spitzer.caltech.edu/spitzermission/observingprograms/legacy/ ** http://ssc.spitzer.caltech.edu/irs/features/ 
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Fig. 2. Top panel - Average spectrum of the equivalent field of view 
of the SH map of the nuclear region of NGC 4945, as seen in the SL 
map. The SL2 order is shown with a dashed line. The most significant 
fine structure emission lines and PAH features in the SL wave range 
are labelled, in particular the ones that overlap with the SH waverange. 
Bottom panel - Average spectrum of the whole field of view of the SH 
map. The higher spectral resolution of the IRS/SH module allows to 
distinguish emission lines like [Ne V] 14.32yum and [CI II] 14.37/jm, oth- 
erwise blended in the IRS/SL module. These spectra are not corrected 
for extinction. 



in Fi sure s l5l l6l [TOl and [TTI are the average fluxes in the 2x2 aper- 
tures. Using the 2x2 pixel extraction box the S/N of the spectra 
is increased by a factor 2 (from the y/A- pixels), and reduces the 
effect of the mismatch between orders mentioned above. This 
corresponds to a standard procedure recommended to obtain IRS 
maps of extended sources. 

The average Full Width Half maximum (FWHM) of the 
IRS/SH Point Spread Function (PSF) has been estimated to 
be ~ 5.16" and ~ 4.31" for the direction along the spec- 
tral dispersion axis and the spatial axis of the slit, respectively, 
while the coiTesponding average PSF FWHM values of the SL2 

1 mostly used in this work) ar e ~ 3.76" and 3.27", respectively 
Pereira-SantaelLa et al.l l2010. their appendix A). The 1-cr devi- 
ation around the average centroid position in SL modules was 
found to be ~ 20 times smaller than the SL pixel size (1.85"). 
While the maximum l-cr deviation in the SH module was found 
to be ~ 0.6", which is considered small compared to the SH 
pixel size (2.26"). 

With the aim to determine whether our line emission maps 
are resolved, we retrieved a partial map of 3x3 SH slit pointings 
of the star P Cygni (AOR 13049088). The map consists of 5" 
steps parallel to the slit and 2" steps perpendicular to the slit, 
similar to the map spacing we have for NGC 4945. These point- 
ings lead to a 10x3 pixel map. In order to test the continuum 




14.0 14.2 14.4 

Rest wavelength \ [firn] 

Fig. 3. Zoom into the [Ne V] complex wavelength range from the av- 
erage IRS/SH spectrum of Fig. |2] At these wavelengths we find the 
first bending mode of HCN at 14.02 yum, and the AGN tracer [Ne V] at 
14.22 pra, almost blended with the fine-structure emission line [CI II] at 
14.37 pm, on the slope of the PAH feature around 14.22 pra. The dotted 
line shows the wavelength at which we found an order mismatch, proba- 
bly producing the unknown feature observed next to the HCN 14.02 pm, 
which could be affected. 



drop off along the (10 pixels) slit, 2x2 aperture spectra were ex- 
tracted in the same way as we do for the SH maps of NGC 4945. 
Then the continuum flux at 15 pm was measured in each posi- 
tion, resulting in a mostly Gaussian variation of the flux along 
the slit. The FWHM ranges from 6 to 7 arcsec between 10 and 
19 microns, far less than the ratio 19/10 expected over this range 
for the PSF of a diffraction limited telescope. The same proce- 
dure was applied to the 12x 6 map of the star KsiDra used by 
iPereira-SantaeUa et all (12010) (AORs 16294912 and 16340224). 
In this case, a FWHM in the range 5.8" - 6.4" was found in both 
the parallel and the perpendicular (to the slit) directions. 

From the SH 2x2 aperture average spectra we obtained maps 
of the starburst tracers [S IV] at 10.51 pm, [Ne II] at 12.81 //m, 
[CI II] at 14.37 pm, [Ne III] at 15.56 pm, [S III] at 18.71 pm, 
as well as the AGN tracer [Ne V] at 14.32 pm, and the pure 
rotational molecular hydrogen lines, H2 S(2) and H2 S(l) at 12.3 
pm and 17.0 pm, respectively. For each line we first estimate 
a continuum level by fitting a second-order polynomial to the 
base of the emission line, or a cubic spline with selected anchor 
points (pivots) when the emission line was on top of a complex 
PAH or silicate feature. We applied the same aperture reduction 
procedure for the SL map, from where we obtained a map of 
the silicate absorption feature at 9.7 pm and the corresponding 
silicate strength as described in Sec. 13.41 We also obtained a map 
of the the molecular hydrogen line H2 S(3) at 9.66 pm. 

3.1. The spectral complex around [Ne V] 14.32pm 

The complexity of the spectral range around the [Ne V] 14.32 pm 
is shown in Fig. [3] An unknown feature at ~13.95 pm seems 
to be present just next to the first bending mode of HCN 
at 14.02 pm. Depending on the temperature of the suiTound- 
ing molecular gas, HCN can be vibrationally excited by ab- 
sorbing the infrared photons at 14.0 pm. This will produce 
a subsequent cascade process that can enhance HCN rota- 
tional lines in the (sub-)millimeter range. This corresponds 
to the IR-pumping scenario proposed to explain the bright 
HCN 7 = 1 — > (and higher) transition observ ed in some 
ULIRGs and Galactic star-forming regions (e.g.. 



11995^ iGarcia-Burillo et al.ll2006t iGuefin et al.ll2007t 



Aalto et al. 



Aalto et al. 
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l2007t iPerez-Beaupuits et alJlIoToh . Although, due to the order 
mismatch (mentioned in Sec.O observed between ~ 13.90 ;um 
and ~ 14/im, also seen by Pereira-Santaella et al. (2010), we are 
unable to conclude whether we detect or not the HCN 14.02 fim 
feature, either in emission or in absorption. 

A larger PAH feature is observed around ~ 14.22 fim, and 
on its rightmost (and less steep) slope the fine-structure emis- 
sion lines of [Ne V] 14.32 fim and [CI II] 14.37 yum are found. 
Chlorine, which has an ionization potential of 12.97 eV, is 
known to play an important role in characterizing the neutral 
gas components in the ISM. When H2 is abundant (optically 
thick in the FUV), it reacts (exothermically by 0.17 eV) with 
[CI III to form HC1+, which leads to the formation of [CI I] and 
H I (|Jurdll974l:[jura & Yorklll978h . This means, chlorine is pre- 
dominantly ionized in H I regions while it is predominan tly neu- 
tral when co ld H2 components are present ( Sonnentruck er et all 
I2OO2I l2003h . The [CI II] 14.37 fim fine-structure emission Une 
was clearly detected (>3cr) in the whole region mapped with the 
IRS/SH module, and presents a similar (although more spread- 
out) distribution than that of the [Ne II] 12.81 fj.m line. 

Because its ionization potential (97.1 eV) is too large to al- 
low production by main-sequence stars, the fine-structure line 
[Ne V] 14.32 fim is commonly used to probe the narrow l ine re- 
gion of AGNs (e.g , lMoorwood et alJll996aHGe nzel et alll998t 
'Sturm et al.' 2002; 'Armus et al .l 120071 : lAlonso-Herrero et alj 
2009; Baumetal. 2010; WillettelalJ |2010l and references 
therein). It has recently been used as a diagnostic tool to un- 
ambiguously identify AGN in galaxies t hat have not been identi- 
fied a s such using optical spectroscopy dGoulding & Alexandeij 
I2009i) . 

3.2. Estimating the iine emission flux 

The emission fluxes were estimated by fitting a Gaussian pro- 
file to the fine-structure lines and integrating the flux above a 
local continuum, which was estimated from a cubic spline inter- 
polation of selected anchor points (pivots). Since we we want to 
determine the actual spatial distribution of the emission lines in 
the nuclear region of NGC 4945, we set a strict 3cr level detec- 
tion for all the fine-structure emission lines, and at least two data 
points (or a full width of ~560 km s 'at 14.32 yum) in the profile 
to consider a feature as a real emission line in the spectrum. For 
a robust estimate of the local continuum we adopted an average 
flux of five points around each of the selected pivots used in the 
cubic spline interpolation. Figure|4]shows the Gaussian fit of the 
[Ne V] 14.32 fim and [CI II] 14.37 fim fine-structure emission 
lines at three different positions of the IRS/SH map. The uncer- 
tainty of each spectral data point is indicated by the error bars. 
The pivots used for the cubic spline interpolation of the local 
continuum (baseline) are shown with filled circles. These piv- 
ots are at the wavelengths 14.201 fim, 14.247 fim, 14.281 fim, 
14.428 yum, 14.531 fim, 14.667 yum, and 14.780 yum. 

The top panel of Fig.|4]shows the spectrum in a region of the 
map with relatively high uncertainties (~1.4% on average) in the 
spectral data and with high RMS (as computed from the base- 
line, and indicated in the legend of the plots) that leads to less 
than a 3cr detection for [Ne V], but a clear detection (>3cr) for 
the [CI II] line. The spectrum of the middle panel was extracted 
from the central 2x2 pixel aperture of the map, which has a S/N 
ratio ~9 in the line and >100 in the continuum, and is among 
the highest S/N in the whole IRS/SH map. This spectrum shows 
very low uncertainties (~0.06% on average) in the data, and clear 
detections of both lines. Although with low uncertainties in the 
data (~0.13%), the bottom panel of Fig. |4] shows less than a 3cr 
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Fig. 4. Gaussian fit of the [Ne V] 14.32 yum and [CI II] 14.37 fim fine- 
structure emission lines at three difi'erent pixels of the IRS/SH map. 
The error bars indicate the uncertainties of the spectral data. The filled 
circles indicates the pivots used for the cubic spline interpolation of 
the baseline. The top panel shows less than 3cr detection for [Ne V] 
but a clear detection of the [CI II] line. The spectrum of the middle 
panel shows clear detections of both lines. The bottom panel shows 
<3cr detection for [Ne V], in spite of the low uncertainties in the data 
points. 



detection for [Ne V] (with a S/N~1.8), but a strong detection of 
[CI II] 14.37 //m. With this Gaussian fit procedure, we produced 
9x9 pixels maps (as described above) of several fine-structure 
lines, PAH features and the silicate absorption feature at 9.7 jum 
(discussed in Sec. 13.41 ). covering a region of about 23" x 23" 
(about 426x426 pc^). The SH maps are shown in Figs. |5] and 
|6] In the case of the [Ne V] 14.32 pm^ the pixels with surface 
brightness lower than 0.1 x 10"'^ W cm"^ sr"' correspond to a 
detection level < 3cr. 

For the particular case of the [Ne V] line, a Gaussian fit of 
the uncorrected for extinction map indicates that this emission is 
only marginally resolved along the major axis, given its FWHM 
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Fig. 5. IRS/SH surface brightness (10 '^ W cm ^ sr ') maps of the fine-structure lines (left panels, from top to bottom) [Ne V] 14.32 yum, 
[Ne III] 15.55 yum, [Ne II] 12.81 jim, (right panels, from top to bottom) [S IV] 10.51 pm, [S III] 18.71 pm, and [CI II] 14.37 pm. The peak fluxes 
are shown with a white pixel. These fluxes are not corrected for extinction, and each of the 9x9 pixels correspond to the fitted flux in the 2x2 
aperture extracted spectrum as described in Sections [2.31 and \T2\ The dashed-line circle represents the ±1" pointing accuracy of Spitzer, and the 
contour lines are labelled. The reference (Aa = 0,A6 = 0) is marked with a cross a nd corresponds to the position R.A.(J2000)=13:05:27.477, 
Dec.(J2000)=-49:28:05.57 of the H2O mega maser reported bv lGreenhill et alJ (Il997l) . 



of 6.4". Along the minor axis the FWHM is 7.3". After extinc- 
tion correction (discussed in Secs. l3.3l and l3.4l i. the [Ne V] emis- 
sion appears unresolved. The other emission lines (uncorrected 
for extinction) are well resolved (FWHM> 7"). 



In order to compare with our data, we produced an image of 
Ay at the same resolution of the IRS/SH maps, using the HST- 
NICMOS H - K colour image. First, we use the same definition 
of foreground screen extinction given in Marconi et al. ( 2000) 



3.3. Visual extinction 

Lower limits on the reddening in the nuclear region of 
NGC 4 945 were obtained fr om the HST-NICMOS H-K color 
image JMarconi et alj|2000l) . An average color H-K - 1.1, 
yielding an Ay ^11 mag, was observed in the region of the 
Paa ring. This Ay is comparable with extinctions larger than 
13 mag, as estimated from the BTa/BiB and PaafHa ratios 
(iMoorwood & 01ivalll988l:lMarconi et aljl2000l) . 



Ay(H - K) 



E{H - K) 
c(H)-c{Ky 



(1) 



where E(H - K) is the colour excess that can be obtained 
from the difference between the observed and intrinsic colours, 
E(H - K) ^ (H - K)- (H - K)q. We use the average intrin- 
sic colour (H - K)o ^ 0.22 ±0.1 mag of spiral and elliptical 
galaxies reported by Hunt et al. (1997). The c{H) and c{K) co- 
efficients are derived from the extinction law, A^ - c{A)Ay. But 
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Fig. 6. IRS/SH surface brightness (10"'~ W cm"~ sr"') maps of the molecular hydrogen lines (left panels, from top to bottom) H2 S(l) and H2 
S(2) at 17.0 yum and 12.3 /jm, respectively, and the H2 S(3) line at 9.7 pm from the IRS/SL map. The right panel shows (from top to bottom) the 
average continuum flux density (10"'^ W cm"- /jm"' sr"') between 14.5 pva and 15.0 yum, and the surface brightness (10"'" W cm"^ sr"') of the 
PAH features at 11.3 yum and 12.7 yum. The details in the maps are the same as in Figure[5] 



in stead of using the A^ ex {Ajl fimY^''^ law {A> I jjin) assumed 
by Marconi et al. (2000), we used the extinction law for the local 
ISM reported by Chiar & Tielen s (2006), which considers solid 
and porous spheres and a continuous distribution of ellipsoids 
in the extinction profiles used for the amorphous silicates in the 
9.7 fim absorption feature. We interpolated the c coefficients at 
A = 1.606 fim and A = 2.218 m n of the HST-NICMOS H and 
K bands, respectively. In lChiar & Tielens (2006) the local ISM 
continuum extinction was found to be described by the expres- 
sion 



logiAJAK) = 0.65 - 2AlogiA) + 134log(Ay 



(2) 



using Ak/A\/ - 0.09 to normalize the extinction to the K band 
(Whittet 2003). Beyond 8 yum the silicate pr ofile of WR-98A 
is sup erimposed using Av/t(9.7 yum) - 18 (JRoche & AitkenI 
119841) . So we used eq.© and Ax/Ay - 0.09 to estimate the 



c(H) and c(K) coefficients for the stellar light-based extinction 
Av{H - K). 



The top panel of Fig.|7]shows the H ST-NICMOS H-K mag 
image reported bv lMarconi et al.1 (l2000l) with the Spitzer-IRS SH 
grid overlaid. We took all the pixels of the H-K image (which 
have higher resolution than our maps) that fall within a 2x2 pixel 
aperture of the SH grid map, and computed the average Ay for 
that particular pixel. After doing this for every single pixel in the 
SH grid we obtain the foreground extinction (at the resolution of 
the SH map) shown in the bottom panel of Fig. |7] The SH pixels 
that cover less than 400 pixels in the H-K image are shown in 
dashed lines. The peak of this stellar light-derived extinction lies 
about 2.3" (one pixel) northeast of the H2O mega maser 
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Fi g. 7. Top panel - HST -NICMOS H-K colour image (mag) reported 
by Marconi et al. (2000). The Spitzer-IRS SH grid is overlaid to show 
the difference in resolution and the different orientation of the fields of 
view. Bottom panel - Foreground screen extinction Ay (-W-A^) estimated 
from the H-K colour image at the resolution of the IRS/SH map. 



3.4. The deep silicate absorption around 9.7 nm 

From the IRS/SL cubes we can study the amorphous siUcate 
grains, which present a broad and deep absorption feature around 
9.7 jum. The presence of PAH emission along the same line of 
sight as the silicate absorption makes it hard to measure the ap- 
parent silicate depth. If one assumes the PAH emission to be 
foreground to the silicate absorption, then a different apparent 
silicate depth is found than when it is assumed that they are 
mixed. Silicate absorption can also be foreground to the PAH 
emissi on. Here we use the method proposed by ISpoon et aP 
(l2007l their Fig. 2) for absorption-dominated spectra to infer the 
apparent strength of this silicate feature. First, we extracted 1-D 
spectra from the SL map rebinned to the 2x2 pixel aperture of 
the SH map described in Sec. [3] For each spectrum, we adopted 
a local mid-infrared continuum by interpolating the feature-free 
continuum pivots at 5.3 jum and 13.55 /vm. 

Because the S/N is not uniform for all the pixels of the SL 
map, and the features around 5.3 //m change depending on the 
proximity to the nucleus of NGC 4945, we used a similar pro- 
cedure as the one used for [Ne V] 14.32 //m (Sec. 13. II ). The flux 
assigned to the pivot at 5.3 //m was adopted as the minimum 
flux observed between 5.1 ^m and 5.5 yum. In order to avoid 
contamination from the H2 S(3) 9.66 /im molecular line, the flux 
estimated for the deepest point adopted at 9.85 /vm was the me- 
dian flux between 9.7 jum and 10.0 /^m. The flux estimated for 
the pivot at 13.55 jt/m was the median flux observed between 
12.9 yum and 14.0 yum. Then we computed the ratio between the 




Fig. 8. Top panel - IRS/SL map of the apparent silicate strength 
5sii estimated at 9.85 pm. Bottom panel - Silicate-based extinc- 
tion Ay (9. 85 /jm) map estimated from the average visual extinc- 
tion t o silicate optical depth ratio Av/t(9.7 fim) = 1 8 for the local 
ISM JRoch e & Aitkenlll984l: iRieke & Lebofskvlll985h and assuming 
t(9.7 pm) a 550(9.85 pm). The contour lines of the stellar light-based 
extinction A\j(H - K) of Fig. [7] are overlaid on the Av(9.85 yum) map 
(colour pixels). 



observed flux (/obs) and the estimated continuum flux (/ct 
9.85 jim, and we obtained the sihcate strength 5 sii as 



5sii = In 



/obs(9.85yum) 



/eont(9.85jt/m) 



,t)at 



(3) 



For sources with a silicate absorption feature, 5 ,11 can be in- 
terpreted as the negative of the apparent silicate optical depth 
(t"9.7 /im)- The top panel of Fig. [8] shows the IRS/SL map (re- 
binned to the 2x2 aperture of the IRS/SH map) of the appar- 
ent silicate strength at 9.85 jum. The siUcate-based extinction 
Av(9.85 pm), shown in the bottom panel of Fig. [8] was es- 
timated from the average visual extinction to silicate optical 
depth rati o Av/t(9.7 pm) = 18, which is a ppropriate for the lo- 
cal ISM ( jRoche & AiflienllT984t iRieke &Lebofskv.l985.) . The 
spatial distribution of the siUcate-based extinction is similar to 
that of the stellar light-based extinction Ay(H - K) estimated in 
Sec l3.3l The peak extinction is found at the same relative po- 
sition (A R.A. K 2", A Dec ^ I "), at about 2.3" (on e pixel) 
northeast of the H2O mega maser (iGreenhill et alJfl997h . 

It is known that extinction estimated from optical or near- 
infrared observations generally underestimates the actual extinc- 
tion if the environment probed is optically thick at the emission 
lines observed. If most of the emitting region is obscured, as 
in the case of the nucleus of NGC 4945, the visual extinction 
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estimate is representative of the surface of the obscured region 
and not the region itself. This effect is reflected in the differ- 
ent extinctions derived from the H - K optical image and the 
silicate-based estimate derived from our mid-IR observations, 
where the peak extinction is about 7 times stronger than that of 
the stellar light-based estimate. Our Av(9.85 /im) is also a fac- 
tor ~1.7 higher than the extinction (Ay = 36;^ jj mag) previously 
inferr ed from ISO observations of the [S III] 18.7/33.5 yum line 
ratio dSpoon et ali2 000). This difference can be explained by the 
peaked nature of the silicate-based extinction map and because 
of the larger apertures (> 14" x 20") of the ISO observations, 
which averages out the extinction to the lower value. 

Figure |9] shows the contour lines of the IRS/SH surface 
brightness map of the molecular hydrogen line H2 S(I) 17.0 yum 
(top left), H2 S(2) 12.3 ^im (top right), the IRS/SL map of the H2 
S(3) 9.7 jum Une (bottom left), and the [Ne II] 12.81 yum line (bot- 
tom right), overlaid on the silicate-based extinction Av(9.85 pm) 
map. The molecular hydrogen emission avoids the obscured nu- 
cleus and peaks between one and three pixels away from the 
highest obscuration. Although the H2 S(2) 12.3 pm (Fig.|6]l has 
a similar distribution as the H2 S(l) line, its peak emission lies 
closer (~ 2.3", one pixel) to the peak obscuration and to the H2O 
mega maser than the other H2 lines. 

Since the distribution of the silicate-based extinction 
Av(9.85 pm) is similar to that of the stellar light-based ex- 
tinction Ay(H - K) we think that the mid-IR derived extinc- 
tion towards the starburst ring may be accurate to correct the 
fluxes of the starburst tracers. That is, the elongated shape of 
the silicate map indicates that the obscuration is mostly asso- 
ciated with the starburst ring rather than with the line of sight 
to the AGN. Hence, the Av(9.85 pm) map can be used to per- 
form an extinction correction for species whose emission emerge 
from the disk. Species associated with the AGN (BLR and 
NLR gas) may be suffering more extinction than the disk does. 
Therefore, the corrected flux derived for the [Ne V] 14.32 pm 
corresponds only to the best lower limit we can derive from the 
current data. Figures [10] and [TT] show the same set of maps as 
in Figs.|5]and|6l but corrected for extinction using the estimated 
Av(9.85 pm) map. Note that after correcting for extinction, all 
the fine-structure emission lines, the average 14.5-15.0 jum con- 
tinuum, and the PAH features peak at about the same position, 
that of the H2O mega maser. The corrected for extinction H2 
lines (Fig. [TT]) show an offset of ~ 2.3" (one pixel) with re- 
spect to the water maser. However, there is a difference of only 
< I" between the centroids obtained from a two-dimensional 
Gaussian profile fit of all the lines. 

The integrated flux densities of the co-added spectra from the 
whole 10x10 aperture of the IRS/SH map (described in Sec.[3]l 
are summarized in Table [T] In this work we include only the 
fluxes of the most prominent emission hnes. These flux es are 
larger than those reported by iBemard-Salas et all (12009 ). be- 
cause our 10x10 aperture is larger than the SH staring aperture. 
In Table[T]we also include the fluxes corrected with the silicate- 
based extinction Av(9.85 pm), and the line widths (FWHM) as 
estimated from the Gaussian fit. 



3.5. Rotation in ttie nuclear region 

Even though the spectral resolution (R~600) of the IRS/SH 
module is relatively low in comparison to the resolution used 
in most kinematic studies, we used the SH spectra to determine 
shifts in the velocity of various lines. Pereira-Santaella et al. 
( I2OIOI) studied the validity of the SH velocity fields by using 
synthetic spectra. They found that the distortion of the wave- 



Table 1. LINE FLUXES" FROM THE lOx 10 PIXEL (-426x426 pc^) 
APERTURE IRS/SH MAP 



Line 



[/im] 



Flux '^ 
[10"-' Wcm- 



Flux(Av)'' FWHM' 

[Ur^'Wcm--] [kms"'] 



[SIV] 



H, 0-0 S(2) 



[Ne II] 



[Ne V] 



[CI II] 



[Ne III] 



H2 0-0S(1) 



[S III] 



10.5105 



12.2786 



12.8135 



14.3217 



14.3678 



15.5551 



17.0348 



18.7129 



1.61 


22.03 


382.68 


0.08 


1. 12 


30.54 


154.07 


889.94 


525.66 


1.99 


11.57 


6.23 


1201.08 


6026.69 


615.54 


118.85 


563.35 


32.37 


5.34 


27.93 


584.70 


0.53 


0.64 


15.44 


12.82 


66.75 


502.97 


0.10 


0.53 


4.04 


127.01 


741.49 


640.81 


2.79 


18.03 


12.74 


252.60 


1781.25 


549.60 


7.64 


63.90 


18.12 


74.87 


618.00 


562.71 


1.37 


10.90 


9.59 



'' These are the fluxes obtained from the co-added spectra of the whole 10x10 aperture 

of the SH map described in Sec.l3] 
^ Rest wavelength of the lines. 
'^ The flux densities of each line are given in the first row. The row below shows the 

conesponding uncertainties. 
'^ Flux density corrected for extinction using an average extinction value Ay - 36 mag 

estimated from the silicate-based Av(9.85//m) (Sec. 13.41 . For [Ne V] this extinction 

correction may be insufficient. See Sec. 14. 41 
^ Line width in km s" obtained from the Gaussian fit. The line profile is not corrected 

for the instrumental profile, which has a FWHM of 500 km/s. 



length scale introduced by the telescope pointing uncertainties 
limits the accuracy of the velocity estimates up to ±10 kms"', 
regardless of the S/N of the spectra, for sources which are neither 
point sources nor uniformly extended sources (e.g., the nuclear 
region of NGC 4945). Considering as well the uncertainty in the 
absolute wavelength calibration (~ 10% of a pixel), which does 
depend on the S/N of individual spectra, iPereira-Santaella et al.l 
estimated uncertainties of 10 - 30 kms"' in the SH velocity 
fields, and concluded that variations of > 20 kms ' in the ve- 
locity maps are likely to be real. 

The derived velocity fields of the [Ne II] 12.81 pm, [Ne III] 
15.56 pm, H2 0-0 S(l) 17.0 pm and S(2) 12.3 pm, and the [CI 
II] 14.37 pm lines are shown in Fig. [12] These correspond to the 
lines with the highest S/N (> 100) in the SH spectra throughout 
the whole region mapped, and therefore the velocity fields de- 
rived from them are considered to be reliable. We also include 
the velocity field of the [Ne V] 14.32 pm line, although with a 
much lower S/N ranging from ~ 1 to ~ 100 throughout the map. 
The pixels with less than a 3cr detection level in the [Ne V] line 
(mostly found at the edges of the map, where the S/N is the low- 
est) are shown as hatched. 

Given that the galaxy major axis runs at PA~ 45°, and the 
fact that the inclination of the nuclear disk has been estimated 
to be smaller (~ 62") than the inclination (~ 80") of the large- 
scale galactic disk dChou et al.ll2007l) . we find that the range of 
velocities probed in a region of ~ 23" x 23" (~ 426 x 426 pc^) 
with the [Ne II], [Ne III], and H2 S(2) lines agree with the range 
of velocities probed by the ISAAC long slit spectra of the H I 
Pf/3 and H2 0-0 S(9) Hnes by Spoon et al. (2003, their Fig.7). 

The velocity structure of [Ne II] (top left panel in Fig. [T2] ) 
shows a gradient along the N-S direction rather than NE-SW, as 
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Fig. 9. IRS/SL map of the silicate-based extinction Av(9.85 yum) with the contour lines of the IRS/SH map of the surface brightness (in units of 
10"'^ W cm"- sr"') of the molecular hydrogen lineH2 S(l) 17.0 yum (top left panel), the IRS/SH map of the H2 S(2) 12.3 ;/m line (top right panel), 
the IRS/SL map of theH2 S(3) 9.7 ;um line (bottom left panel), and the IRS/SH map of the [Nell] 12.81 pm (bottom right panel). The figures show 
that the H2 emission is stronger in a region >2.3" away from the peak obscuration, although the H2 S(2) is the closest to the peak obscuration and 
to the H2O mega maser, as well as the [Ne II], which peaks at the H2O mega maser (within ~ 1"), as seen in Fig. [5] 



in the case of H2 0-0 S(l) (top right panel in Fig.[T2l). A simi- 
lar difference in velocity structure is observed in the [Ne III] and 
H2 0-0 S(2) velocity maps (middle left and middle right panels 
in Fig. [12] respectively). These different kinematic components 
may be caused by different origins related to the species emit- 
ting, or by different levels of extinction (extinction law at differ- 
ent wavelengths), which allow us to probe deeper regions at one 
wavelength than another. This is further discussed in Sect. 14.21 



On the other hand, the relatively uniform central region of 
the [Ne V] velocity field may be due to the lack of spectral reso- 
lution and the lower (factor 1 to 100) S/N level obtained in this 
line. This implies that we are unable to spectrally resolve, at a re- 
liable level, the rather small (< 50 km s"') velocity shifts in the 
[Ne V] line. Nevertheless, with a higher spectral resolution we 
would expect a weaker rotation of [Ne V] in comparison to the 
rotation shown by the starburst tracers (from about -120 km s"' 
to ~ 100 km s"'). If the [Ne V] is exclusively related to the AGN 
NLR, this would mean that the NLR gas does not feel the grav- 
itational pull of the large gas mass interior to the starburst ring. 
But it would feel the pull from the SMBH, showing not neces- 
sarily the same kinematics as the gas in the starburst ring. In fact, 
the group of three pixels next (~ 43 pc) to the adopted location 
of the AGN (the H2O mega maser), shows positive velocities (in 
the range ~ 10-70 kms"' with uncertainties of <1% in the 
line fit) which is opposite to the characteristics of the velocity 
fields of the starburst tracers. This might be a sign of the kine- 
matically decoupled component discovered at the center of the 
disk with interferometric maps of the 7 = 2 — > 1 transitions 



of I^CO, i^CO, and C'^O (" Chou et alJl2007l) . However, higher 
spectral and spatial resolution observations of the [Ne V] line 
are required to reliably conclude on this. 



4. Discussion 

4. 1 . Tracing the circumnuclear starburst ring 

The [Ne II] and [S III] fine-structure lines are our cleanest tracers 
of star formation (H II regions). If an AGN is present, however, 
the [Ne III] and [S IV] lines will have contributions from both the 
starburst and the AGN. In NGC 4945 the AGN contamination to 
[Ne III] is likely low, given the faintness of [Ne V]. In fact, at 
the peak flux of the [Ne V] and [Ne III] maps (Fig. [5]) the former 
corresponds to only about 5% of the [Ne III] flux. 

From Table [T] we derive the respective luminosities of 
L[N,v] =7.42xl02''Wsr-i andL^Neiii] = 1. 76 x 10^' W sr"', us- 
ing a luminosity distance Dl-3.^2 Mpc. However, from the re- 
lation between L[NeV] and L[Nem] bv iGoriian et al.l (l2007h (which 
holds for a group of Seyfert 1 and 2, AGNs, ULIRGs and ra- 
dio galaxies), the predicted [Ne III] luminosity should be ~ 
3.28 X lO"*" W sr"', that is a factor -5.4 lower than observed. 
This means that NGC 4945 has a slighdy lower AGN contri- 
bution to the [Ne III] line than the outlier NGC 3079 found in 
the sample by Goijian et al. (2007), for which the [Ne III] lumi- 
nosity is a factor ~5 brighter than expected based on its [Ne V] 
luminosity. 
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Fig. 10. IRS/SH surface brightness (10 '^ Wcm ^ sr ') maps, corrected for extinction using the extinction law for the local ISM from 
IChiar & TielensI ( 12006.) and the silicate-based extinction Av(9.85 jim) (Figs [8] and |9}- The left panels show (from top to bottom) the fine- structure 
lines [Ne V] 14.32 ^lm, [Ne III] 15.55 fim, and [Ne II] 12.81 yum. The right panels show (from top to bottom) the [S IV] 10.51 yum, [S III] 18.71 pm, 
and [CI II] 14.37 pm. lines. The extinction corrected fluxes and notations are as in Fig.[5] Note that after correcting for extinction, all these emission 
lines (and the average continuum) peak at about the same position, that of the H2O mega maser. 



The [Ne III] 15.56/[Ne II] 12.81 pm line ratio (Fig. [T3]l 
ranges between ~0.13 and ~0.27 above and below the major 
axis of the region mapped. The difference with the lowest ra- 
tios (0.06-0.13), which are found along the northeast-southwest 
axis, is more pronounced than in the [Ne V] 14.32/[Ne II] 12.81 
pva line ratios {top left panel in Fig.[T3ll. 

The lowest [Ne III] 15.55 jum/[Ne II] 12.81 pm line ratios 
coincide with the 100 pc-scale circumnuclear starburst ring of 
NGC 4945 (bottom panels in Fig. [T4 ]). detect ed in Pao- emis- 
sion with HST NICMOS ( Marconi ^i]|2003). The increasing 
gradient above the starburst ring seems to mim ic the conical- 
shape d cavity traced by the H2 1-0 S(l) line (iJVIarconi et alJ 
l2000h . However, the ratio also increases toward the south-east 
direction from the starburst ring, which rather implies that the 
ring has different properties than its surroundings. 



In fact, a similar increase of [Ne III]/[Ne II] ratios 
awa y from the nucleus is o bserved in the nuclear maps by 
Per eira-Santaella et al.l (1201 Ol their Fig. 5), with exception of 
NGC 7130, the only galaxy classifi ed as LINER/ Seyfert 2 in 
their sample. According to (Sniide rs et alJ 120071) the higher 
the metallicity, the higher the age of the stellar population 
and the higher the gas density, the lower should be the 
[Nelll] 15.55 //m/[Ne II] 12.81 //m ratio. 

The [S III] 18.71/[Ne II] 12.81 Une ratio is considered a good 
tracer of densities in the range 10^ and 10^ cm""* because their 
similar excitation potentials (21.6 eV and 23.3 eV, for [Ne II] 
and [S III], respectively) and significantly different critical den- 
sities (6.1 X 10^ cm"^ for [Ne II] and 1.0 x lO"* cm"^ for [S III]) 
make this ratio less sensitive to the hardness of the radiation 
field than to the density of the ISM. According to the models 
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Fig. 11. IRS/SH surface brightness (10"'^ Wcm"~ sr"') maps, corrected for extinction using the extinction law for the local ISM from 
[Chiar & T ielens (2006) and the silicate-based extinction Av(9.85 pm) (Figs[8]and|9}- The molecular hydrogen lines H2 S(l) 17.0 pm, H2 S(2) 
12.3 yum, and H2 S(3) 9.7 yum (from the IRS/SL map) are shown in the left panel. The right panel shows (from top to bottom) the average contin- 



uum flux density (10 '^ W cm ^ pm ' 



sr ) between 14.5 pm and 15.0 /im, and the the surface brightness (10 W cm - sr ) of the PAH features 



at 11.3 prs\ and 12.7 jum. The extinction corrected fluxes and notations are as in Fig.|5] After correcting for extinction, the [CI II] line and the PAH 
features also peak at about the same position of the H2O mega maser. 



by ISniiders et aP (l2007h a lower [S III]/[Ne II] ratio would in- 
dicate larger densities. In Fig.[T3]the [S III] 18.71/[Ne II] 12.81 
line ratios are also lower in the center than around it. This ra- 
tio ranges between 0.002 and 0.013 in the center, which indicate 
densities larger than 10^ cm""^ in a >5 Myr old starburst sys- 
tem with solar metallicity and relatively high (q = 8 x 10**) 
ionization parameter (Snii ders et al.l 120071 their Fig. 5). Based 
on our observed [Ne III]/[N e II] ratios, the [S III]/rNe II] ratios 
predicted with the model bv lPereira-Santaella et al.l (1201 Ol their 
equation 2) are about ten times larger than the observed ratios. 
These can be a consequence of the about 10 times higher extinc- 
tion found in NGC 4945 than in the sample of galaxies used by 
iPereira-Santaella et all (l2010l) . 



The [Ne III]/[Ne II] line ratios obtained with the extinction 
correction {middle right panel in Fig. [T3T l are just ~9% larger 
than without correction. This relatively small change after the 
extinction correction is because even in a high-extinction situ- 
ation the differential extinction between [Ne III] and [Ne II] is 
small, given that both lines are closely spaced in wavelength and 
not in one of the silicate absorption features. On the other hand, 
because of their larger differential (wavelength) extinction, the 
[S III]/[Ne II] ratios corrected for extinction do change signifi- 
cantly from a factor ~50 in the center (where the extinction is 
larger) to a factor 3 away from the center (where the extinction 
is lower). 

Note that only the ratios from our co-added fluxes (Table [TJ 
can be compared to other galactic nuclei, as the 10x10 aper- 
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Fig. 12. The left panel shows (from top to bottom) the velocity fields of the fine-structure lines [Ne II] 12.81 jiva, [Ne III] 15.56 jim, and 
[Ne V] 14.32 /jm. The right panel shows (from top to bottom) the molecular hydrogen lines H2 0-0 S(l) 17.0 yjm, and S(2) 12.3 yum, and the [CI 
II] 14.37 yum. The maps, covering a region of ~ 23" x 23" (~ 426 x 426 pc^), are consistent with a rotating disk, with the exception of the [Ne V] 
line, which present a relatively uniform (most likely not spectrally resolved) velocity field. 



ture is comparable to the size probed in any of the more distant 
galaxy nuclei. The [Ne III]/[Ne II] ratio at the position of the H2O 
mega maser is about 10% lower than the ratio obtained from the 
fluxes of the 10x10 co-added spectrum (Table[TJ. 

Since in starburst environments the [Ne II] and [Ne III] 
emission lines are expected to be driven mainly by photo- 
ionization (e.g. jHo & Ketolf2007h . the lower ratios found along 
the northeast-southwest axis are likely due to a [Ne II] emission 
enhanced by the starburst ring. Although, these low ratios are 
also cons istent with a ratio [Ne III]/[Ne II]<0.1 found in shocks 
dBinette et al. 1985) where t he low ionization line [Ne II] can 
also be enhanced (I\^ll992h . 

On the other hand, the highest [Ne III]/[Ne II] ratios found 
above and below the starburst ring are larger than those typ- 
ically found in shocks. These higher ratios could actually 
be tracing an additional contribution to the [Nelll] emission 



by a conically shaped narrow line region (NLR). Previously, 
no evidence was found for the existence of such a NLR, 
given the abs ence of 5007 A [O III] emission in the central 



AK 



800pcx800pc dMoorwood et alJl996bl) . Instead the conical cav- 
ity trace d by optical a nd near-inf rared Une and con tinuum trac- 
ers ( Mo orwood et al J [l 996b: Ma rconi et al.l 120001) was associ- 
ated with a starburst super wind. However, given the high extinc- 
tion implied by our observations, optical [O III] emission can be 
easily attenuated; mid-infrared [Ne III] emission less so. A hypo- 
thetical conical NLR would be difficult to identify from a [Ne III] 
map alone due to blending with the nuclear starburst component. 
But it may manifest itself more clearly in a [Ne III]/[Ne II] map, 
in which the contribution of the starburst component diminishes 
quickly beyond the nucleus proper. 

The [Ne III]/[Ne II] -0.06-0. 13 ratios observed along the 
starburst ring are in close agreement with the ratios observed 
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Fig. 13. Left panels - Integrated flux ratios between the fine-structure lines (from top to bottom) [Ne V] 14.32/[Ne II] 12.81, 
[Ne III] 15.56/[Ne II] 12.81, [S III] 18.71/[Ne II] 12.81, and [S IV] 10.51/[S III] 18.71. Right panels - Same as above, but considering the 
fluxes corrected for extinction. The hatched areas correspond to pixels with a <3cr detection level in the faintest emission lines [Ne V] 14.32 /im 
and [S IV] 10.51 //m. 



in some galaxies of the ISO starburst sample and are consistent starburst of 5 - 8 Myr jThornlev et al.ll200u , their Fig. 6). In the 
with burst timescale models predicting a relatively short-lived whole region mapped, we observe ratios [Ne III]/[Ne II]<0.3, 
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Fig. 14. HST NICMOS i mage of Paor (in uni ts of IQ--' W cm-^) of 
the nucleus of NGC 4945 dMarconi et al.ll200(]h . The not corrected for 
extinction contours (labelled in units of 10"'^ W cm"^ sr"') are the fine- 
structure lines [Ne II] 12.81 /jm (top panel) and [Ne V] 14.32 pm (mid- 
dle panel). The [Ne V] flux peaks at about one pixel to the north- 
west of the H2O mega maser (Fig. |5]l, The bottom panel shows the 
[Ne III] 15.55/[Ne II] 12.81 /jm line ratio not corrected for extinction. 
The ratio is lower along the starburst ring traced by the Paa line, while 
it increases away from the center 



which are a factor ~3 lower than the ratios observed in the 
sample of quasars reported by Veilleux et al. (2009). However, 
the ratio [Ne III]/[Ne II]~0.07 observed at the nucleus of 
NGC 4945 is comparable to the ratios observed in the nucleus 
of 5 (out of 16) galaxies in the sample of LIRGs studied in 
iPereira-Santaella et all (12010 . their Fig. 14). 

On the other hand, the differential extinction between [S III] 
and [S IV] is larger than between [Ne II] and [Ne III]. Because 
[S IV] is sitting close to the deepest point of the silicate absorp- 



tion feature at 9.7 jj.m (Fig©, the [S IV] 10.51/[S III] 18.71 j^m 
line ratios corrected for extinction (bottom right in Fig.[T3Ti are a 
factor ~2 larger than those obtained without correction (bottom 
left in Fig. [13). In both cases the highest ratio is found about 
2.3" (one pixel) northwest of the H2O mega maser. The un- 
coiTected for extinction map-averaged ratios (from the fluxes in 
Table[B of [Ne III]/[Ne II]~0.11 and [S IV]/[S III] -0.02 are simi- 
lar (within 10%) to the ratios obtained from the SH staring obser- 
vations bv lBernard-Salas et all ( |2009|) . They position NGC 4945 
among the sources with the lowest hardness of the radiation 
field in their samp le of starburst galaxies, which, according to 
|Bernard-Salas et al., may be an indication of an old (or small) 
population of massive stars in NGC 4945. However, this can be 
concluded only for the outermost surface that we can probe with 
[S IV] and [S III] lines, and we cannot rely much on the extinc- 
tion corrected data for this analysis, since the [S IV]/[S III] ratio 
depends more on the extinction law used than the [Ne III]/[Ne II] 
ratio. 



4.2. Molecular hydrogen gas 

As can be seen in Figures |5] & IB the emission centroids of the 
H2 0-0 S(l), S(2) and S(3) lines are clearly off'set from the cen- 
troids of the low ionization fine-structure lines, which all peak 
close to the H2O mega maser Such an offset is hard to under- 
stand if the H2 0-0 S(l), S(2) and S(3) lines and the HII region 
gas would be exposed to similar amounts of extinction. The dif- 
ferential extinction between H2 S(2) and [Ne II], at 12.28 and 
12.81 jum, for instance, is very small. The same is true for the 
wavelengths of H2 S(l) and [S III] at 17.0 and 18.71 yum, re- 
spectively. The lack of H2 0-0 S(l), S(2) and S(3) emission at 
the position of the mega maser can be explained if the molec- 
ular hydrogen emission in the nucleus would be more strongly 
attenuated than the HII region gas. That is, the inner H2 emis- 
sion would suffer additional/differential extinction that is not ac- 
counted for by the silicate feature. On the other hand, analysis 
of the higher pure rotational lines observed in the SL module 
(Spoon et al. in preparation) and with VLT-ISAAC ( Spo on et al.l 
[2003, their Figure 7) indicates that the S(5), S(7) and'S(9) fines 
originate progressively closer to the nucleus. This seems hard 
to reconcile with the above scenario, given that A(A)/Av in the 
S(5)-S(9) wavelength range is likely higher than in the S(l)-S(2) 
wavelength range. 

We instead favour a scenario in which the lower pure ro- 
tational lines S(l) and S(2) mainly originate from an unob- 
scured extra-nuclear component associated with the super-wind 
cone as seen in the HST NICMO S map of the vibrational H2 
1-0 S(l) emission jMarconi et al . 2000). Figure [15] shows the 
IRS/SH maps of the pure rotational H2 S(l) and S(2) lines over- 
laid on the H2 1-0 S(l) map. Interestingly, the two pure rota- 
tional lines peak further from the apex of the cone than the 1- 
S(l) line does. The 0-0 S(l) line (upper level energy 1015 
K) more so than the 0-0 S(2) fine (r„p=1682 K). In this sce- 
nario the higher pure rotational lines S(5), S(7) and S(9) would 
originate mainly from the starburst ring or within. There one 
likely finds the high critical densities (nciit > 10^ cm"^) as well 
as high gas temperatures (Tk > 10^ K) needed to thermalize 
these lines (as se en in a ctive gala xies like NGC 1068, Mrk 231, 
and A rp 220; e.g. iPerez-Beaupuits et al.l2007ll2009l:[AaIto et al.l 



2007t lvan der Werf et al.ll2010[) . but fliev wiU also be aifected by 
strong extinction. The latter may cause the nuclear component 
of the 2.12 fim H2 1-0 S(l) fine (r„p=6950 K) to suffer far more 
extinction than the longer wavelength H2 0-0 S(5)-S(9) lines and 
therefore to be associated solely with tfie conical cavity. In all of 
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Fig. 15. IRS/SH contours of the rota tional H; 0-0 S(l) (l eft) and H2 0-0 S(2) (right) lines overlaid on the HST NICMOS map of the vibrational 
H2 1-0 S(l) molecular hydrogen line JMarconi et alj|2000h . 



this the H2 0-0 S(3) line, originating at an upper level energy of 
2500 K, may have both strong nuclear and cavity components. 
Due to the extreme extinction the nuclear component would be 
suppressed, leaving only the cavity component for us to see. 

4.3. Excitation of H2 

For LTE conditio ns, and an ortho-to-para ratio of 3 
(iNeufeld et al.ll20() o. their Fig. 13), we can estimate the excita- 
tion temperature of the molecular hydrogen throughout the re- 
gion mapped with the IRS/SH module, from the ratio between 
the flux of the H2 0-0 S(2) 12.3 yum and H2 0-0 S(l) 17.0 yum 
lines as 



T - -- 



.^S(2) _ .^S(l) 
^ up ^ up 



log(Fs(2)yS(l)As(\)gS(l)) - logiF5(i)Vs(2)As(2)gS(2)) 



K,(4) 



where Fs(i) and Fs(2) are the integrated flux densities of the H2 
0-0 S(l) and S(2) lines, v is the corresponding rest frequency 
(~ 1.75 X 10'^ Hz for S(l) and ~ 2.44 x lO'^ Hz for S(2)), and 
A and g are the respective Einstein A-coefficients and statistical 
weights of each transition. The upper level energy Tup of the 
transitions is in units of K. 

Figure[T6] shows the IRS/SH maps (not corrected for extinc- 
tion) of the H2 0-0 S(2) to S(l) total flux ratio (top left panel), 
and the coiTesponding map of the excitation temperature T^,^ 
(K) (bottom left panel) of the molecular hydrogen, estimated 
from eq.(|4]i. The peak excitation temperature of ~ 528 ± 31 K 
is reached >2.3" (one pixel) to the north of the H2O mega 
maser. The right panels of Fig. [16] show the H2 0-0 S(2)/S(l) 
ratio and the estimated excitation temperature from the respec- 
tive H2 fluxes corrected for extinction with the silicate-based 
Av(9.85 yum). The correction for extinction leads to a peak tem- 
perature of 421 K, which is about 100 K lower than the temper- 
ature derived from the non-corrected fluxes. However, the distri- 
bution of the temperature, and the position of its peak value, do 
not change. 

Since galaxies are complex systems, the H2 gas is not ex- 
pected to be at a single temperature. Besides, a long line of sight 
can probe different excitation environments (e.g., PDRs, XDRs, 
shocks) like in the sample of galaxies studied by Rouss elet alJ 
(l2007h . Previous estimates of the excitation temperature, based 
on SWS observations of the (0-0) S(0) and S(l) fluxes detected 
in the nucleus of NGC 4945, led to a cooler H2 component with 
Tex ~ 160 K (corresponding to about 9% of the total H2 mass), 
while a temperature T^x ~ 380 K (about 0.4% of H2 mass) was 



estimated from the S(l) and S(2) fluxes (S poon et al.ll2000l their 
Table 6), which is in close agreement with our result consider- 
ing the size of the big SWS slit. This excitation temperature is 
similar to the temperature T^x = 36 5 ± 50 K derived fo r the 
low-energy transitions in NGC 1377 dRousselet al.ll2006h . but 
higher than the temperature of T^x - 292 + 6 K estimated from 
the S(0)-S(3) lines observed in NGC 6240 (Higdon et al. 2006i) . 
Since we cannot use the map of the H2 S(3) 9.7 /im line due 
to the flux mismatch between the SH and SL modules mentioned 
in Sec. |3] we would only be able to base our estimate of the 
warm H2 mass on the S(l) and S(2) lines. In fact, even if we had 
been able to use the S(3) line, most of the mass is revealed by 
including the S(0) line. This line has however not been mapped, 
since the spatial resolution of the LH module is much lower than 
that of the SH and SL modules. 



4.4. AGN dominated [Ne V] emission? 

The detection of [Ne V] emission from the nucleus of NGC4945 
is not surprising. The galaxy hosts a 1.4x10^ M0 super 
massive black hole, as revealed by H2O maser observations 
dGreenhill et al.lll9 97) with an intrinsic 0. 1-200 keV luminosity 
of 1.8x10^3 erg s^^ (Gu ainazzi et al. 2000) . 

More remarkable is the faintness of the [Ne V] emission. 
Only 10 '' of the infrared luminosity of the galaxy is detected in 
the [Ne V] line (Table [T]i- This is surprisingly little for an AGN 
which can easily account for the entire bolometric emission of 
the galaxy dMarconi et al.ll2000l) . It is also 100 times less than 
L([Ne V])/LiR for galaxies whose IR lumino sity is dominated by 
AGN activity dGoulding & Alexandeiil2009h . 

Interesting is also the size of the [Ne V] emitting region, 
which can be reasonably well constrained given that at a dis- 
tance of 3.82 Mpc one arcsecond corresponds to 18.5 pc. Since 
the spatial profile of the [Ne V] map (Figure|5]) is only marginally 
wider than the PSF, the radius of the [Ne V] emitting region will 
be less than 55 pc (~3")- This upper limit is large in comparison 
to the size of coronal line regions (e.g. ~19 pc for C ircinus as 
probed by the 2.48/im [Si IV] line: lPrieto et al.ll2004 . but small 
compared to the kpc scale of narrow* hne r egions as probed with 
[O III] (e.g. Circinus: [Marconi et alJll994 . 

Both the faintness and the observed extent of the [Ne V] 
emission are likely affected by strong extinction within the nu- 
cleus. However, if the narrow line region were to extend above 
the plane of the molecular disk in which the AGN and the 
circumnuclear starburst are embedded, the amount of extinc- 
tion would be less than probed by the silicate absorption fea- 
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Fig. 16. Left panels - IRS/SH maps of the H2 0-0 S(2) to H2 0-0 S(l) total flux ratio {top left panel), and the map of the estimated r„ (K) of the 
molecular hydrogen. The peak excitation temperature of ~528 K is reached about >2.3" (one pixel) to the north of the H2O mega maser. Right 
panels - Same as above, but using the extinction corrected fluxes of the H2 S(2) and S(l) lines. The extinction correction leads to an ~100 K lower 
peak temperature. 



ture (Av=50-60). The amount will, however, certainly be less 
than the line of sight extinction to the AG N broad line region, 
Nh=2.4x10^'* cm"^, (iGuainazzi et al.l20doh . which is equivalent 
toAv=1300. 

The extinction probed by the silicate absorption feature is 
consisten t with the F(rNe Vl) /F(14-195keV) ratio found for Syl 
galaxies dWeaver et al.ll2010) when multiplied by the absorption 
corrected 14-195keV flux for NGC 4945. The latter depends 
on the geometry of the absorber, and ranges between 5 and 10 
times the Swift- BAT flux of NGC 494 5 (F(14-195keV)=33±0.8 
xlO"'^ W cm"^: iTueller et all (|2010|) ) for a spherical and edge- 
on disk geometry, respectively (K. Iwasawa, private communi- 
cation). Taking further into account the factor 2.5 scatter in the 
F([Ne V])/F(14-195keV) ratio for Syl galaxies, the attenuation 
on the [Ne V] flux ranges between 12 and 157, which corre- 
sponds t o Av=55-112 m ag for our adopted local ISM extinction 
law of Chiar & TielensI 0200 6). This result would indicate that 
the NLR is buried along with the other components of the nu- 
clear molecular disk. 

The strong obscuration to the NLR in NGC 4945 may serve 
as a warning that mid-infrared NLR tracers such as [Ne V] and 
[O IV] may not always be suitable as tracers of the A GN lumi- 
nosity as suggested bv iGoulding & Alexander! ( l2009h . Galaxies 
in this risk group may be recognized by the presence of a deep 
silicate absorption feature in their mid-infrared spectra and a 
classification 3 A, 3B or 3C in the mid-infrared galaxy classifica- 
tion scheme oflS poon et"an (2007). Besides NGC 4945 galaxies 
in this group include LIRGs such as NGC 4418 and ULIRGs like 
Arp220, IRAS 08572+3915 and IRAS F00183-71 11. 



The example of NGC 4945 may also serve as a warning 
that the non-detection of NLR tracers in AGN does not neces- 
sarily imply the absence of a NLR. Instead it may indicate that 
the NLR is so strongly attenuated that the line does not emerge 
above the continuum noise. In these cases, only deep mid-IR 
observations, like the ones presented here, would reveal these 
tracers. 

In NGC4945 the [Ne V] line was detected at 5% of the con- 
tinuum flux in a spectrum with continuum S/N of 140 in the 
14-15jum range. Other galaxies which may host an AGN, and 
which have been observed at similar S/N (among them ULIRGs 
like Arp220 and IRAS 08572+3915 ), do not show evidence for 
[Ne V] emission (lArmus et al.ll2007h . The NLR in these galaxies 
may hence be even more strongly obscured than in NGC 4945, 
or the AGN may contribute a far smaller fraction of the bolomet- 
ric power than in NGC 4945. The NLR may also simply not exist 
if the ionizing photons needed to form the NLR are 4;r obscured 
close to their origin. 

4.5. [Ne V] emission from supernovae? the template 
Cassiopeia A 

Although the [NeV] emission is unresolved along the major axis 
and only marginally resolved along the minor axis (as discussed 
in Sec.|3]l, the fact that the [Ne V] line has also been observed in 
supernova remnants (SNRs) (e.g. lOliva et al.|[T999t ISmith et al.l 
2009) raises valid questions: could the [Ne V] emission observed 
in NGC 4945 be powered purely by SNRs?, and if so, how many 
SNRs would be needed to reproduce its [Ne V] flux? Below we 
explore the alternative of [Ne V] emission being powered by the 
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SNRs in the nuclear region of NGC 4945. We note, though, that 
the number of supernovae needed depends not only on the total 
[Ne V] flux but also on the actual size of the SNR. 

Our analysis is based on previous Spitzer/IRS SL module ob- 
servations of Cassiopeia A (Cas A) from which the global dis- 
tribution of fine-struc ture lines covering a 5.3x5.3 arcminute^ 
area was presented bv lSmith et all (l2009l) . The total [Ne V] flux 
of Cas A was not reported before due to the difficulties to ex- 
tract the flux with enough reliability. Using fitting techniques 
tailored to the [Ne V] line, we now are able to state its total flux 
as 1.5 + 0.7 X 10"'^ W cm~'^. Conside ring an average extinction 
Av ~ 5 mag dHurford & Fesenlll996h the corrected flux usin g 
the extinction law for the local ISM bv lChiar & TielensI (l2006h . 
is- 1.9±0.9xl0-"* Wcm-2. 

On the other hand, the total [Ne V] flux in NGC 4945 is 
~ 5.3 ± 0.5 X 10"^' W cm"^. In order to estimate the dered- 
dened [Ne V] emission observed in NGC 4945, we computed an 
average extinction value Ay ~ 36 mag from the Av(9.85 fim) 
map (Fig.[8]l. Note that due to the stronger extinction at the cen- 
tral pixels, this represents just a lower limit to the actual extinc- 
tion on the [Ne V] line. Thus the lower limit for the extinction- 
corrected total [Ne V] flux is -27.9x10"^' W cm"^ 

Consider ing a canonical distance of 3.4 kpc for Cas A 
(ISmith et alJ l2009. and references therein) and ~3.82 Mpc for 
NGC 4945, we found that the ratio between the total [Ne V] flux 
of NGC 4945 and Cas A, corrected by the square of the ratio 
between the distances (£>NGC4945/£'casA)^, is about 19x10^. This 
represents a lower limit to the number of contemporaneous Cas 
A type supernovae needed to produce the extinction-corrected 
[Ne V] emission observed in NGC 4945. Another extreme lower 
limit can be obtained assuming no extinction in the NGC 4945 
flux, from which the same exercise indicates that about 4x10^ 
Cas A type SNRs would be needed^ 

According to Lenc & Tingayj (|2009|) . the upper limit to the 
supernova rate in NGC 4945 is 15.3 yr~\ which means the indi- 
vidual supernovae, for the not corrected and corrected for extinc- 
tion lower limits estimated above, need to emit at Cas A levels 
for about 260 and 1240 years, respectively. These are reasonable 
nu mbers given that th e age of Cas A is estimated to be about 330 
yr jFesen et al.ll2006h . while the kinematic age of another [Ne V] 
emitting SNR, E0102, has been estimated to be 2050 + 600 yr 
(iFinkelstein et al. 2006). This implies that the number of SNRs 
estimated above can indeed persist long enough to power the 
[Ne V] emission observed in NGC 4945, at least from the point 
of view of the SNR p roperties used in our estimate. However, 
iLenc & Tingavl (l2009l their section 4.5) also determine a me- 
dian SN age of 85 yr, which translates into a median supernova 
rate of -0.12 SNe/yr, or a much (15.3/0.12-128 times) longer 
radiative phase (-3-16x10"* years) needed for the individual su- 
pernovae. Therefore, the [Ne V] flux observed in NGC 4945 is 
likely not powered by a population of supernova remnants. 



5. Final remarks 

We have mapped the central region of NGC 4945 with the SH 
and SL modules of the Spitzer InfraRed Spectrograph. From the 
SH spectral cubes we produced maps of fine-structure emission 
lines [S IV] at 10.51 A^m, [Ne II] at 12.81 jum, [CI II] at 14.37 ^m, 
[Ne III] at 15.56 fim, [S III] at 18.71 fim, the AGN naiTow-Hne 
region tracer [Ne V] at 14.32 /vm, and the molecular hydrogen 
lines, H2 S(2) and H2 S(l) at 12.3 /urn and 17.0 jum, respectively. 
From the SL spectral cubes we obtained maps of the H2 S(3), 
and the silicate absorption feature at 9.7 fim. 



We present the first map of [Ne V] 14.32 /jxn towards the 
nucleus of NGC 4945 with flux detection levels down to 0. 1 x 
10"'^ W cm"^ sr"' per pixel. We produced and estimated an ex- 
tinction map Av(9.85 fim) based on the apparent strength of the 
9.7 fim silicate absorption feature, which traces the contours of 
the starburst ring at a ~5" spatial resolution. 

Most of the emission lines are found to peak on the nucleus, 
within the uncertainty of the astrometry. Only the warm molecu- 
lar hydrogen emission shows a maximum about 60-100 pc NW 
of the nucleus. After correction for extinction the distribution 
of the H2 rotational emission is more concentrated in the nu- 
clear region, but its peak emission is still slightly offset from the 
peak of the other emission lines, that of the H2O mega maser 
Thus, we favour a scenario in which the lower pure rotational 
lines S(l) and S(2) mainly originate from an unobscured extra- 
nuclear component associated with the super-wind cone as seen 
in the HST NICMOS map of the vibrational H2 1-0 S(l) emis- 
sion, with an intrinsic excitation trend toward the nucleus which 
is reflected in the higher level S(5)-S(9) lines. 

We found that the map-integrated [Ne V]/[Ne II] ratio is con- 
sistent with ratios observed in starbursts rather than in AGNs. 
The [Ne III]/[Ne II] <0. 13 ratios observed along the starburst ring 
are likely due to an excess [Ne II] emission driven by the star- 
burst ring, or to high density (> 10^ cm"^) ISM gas in the cir- 
cumnuclear disk. 

A range of extinction Ay ~ 55-112 mag (which corre- 
sponds to an attenuation of a factor 12-160) estimated for [Ne V] 
from our observed [Ne V] flux and the absorption-corrected 14- 
195 keV Swift-BAT flux, indicate that mid-infrared NLR tracers 
such as [Ne V] and [O IV] may not be trusted as tracers of the 
AGN luminosity of galaxies with a deep silicate absorption fea- 
ture. 

A new estimate of the total [Ne V] flux in Cassiopeia A in- 
dicates that at least 4-19x10"' Cas-A type supernova remnants, 
with ages between >260 (no extinction correction) and >1240 
(on extinction correction) years, would be needed to power the 
[Ne V] emission observed in NGC 4945. However, given the ac- 
tual median age of SNRs observed in NGC 4945, and the uncer- 
tainty in the true extinction of its [Ne V] emission, SNRs are not 
likely to fully reproduce the [Ne V] flux observed in the nucleus 
of this galaxy. 
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